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ABSTRACT 

We present the initial results of a Hubble Space Telescope/ACS snapshot survey of 16 nearby, edge- 
on, late-type galaxies covering a range in distance from 2 to 19 Mpc. The images of these galaxies show 
significant resolved stellar populations. We derive F606W and F814W photometry for >1.2 million 
stars, and present color-magnitude diagrams which show a mixture of young, intermediate, and old 
stars in each galaxy. In one of the fields we serendipitously detect stars from the Large Magellanic 
Cloud. We also identify a candidate young dwarf galaxy lying ~2 kpc above the plane of NGC 4631. 
For the nearest six galaxies, we derive tip of the red giant branch distances and demonstrate that 
these galaxies fall on the K-band TuUy-Fisher relation established in clusters. From the color of the 
red giant branch, we also find evidence that these galaxies possess a metal-poor thick disk or halo 
population. 

Subject headings: galaxies: distances and redshifts - galaxies: dwarf - galaxies: Magellanic Clouds - 
galaxies: spiral ~ galaxies: stellar content - techniques: photometric 



1. INTRODUCTION 

This paper is the first in a series presenting results 
from a Hubble Space Telescope (HST) snapshot program 
of nearby edge-on galaxies. These galaxies are mostly 
low mass, late-type spirals with circular velocities be- 
low 200 km/sec. Although the project was designed to 
study the dust content of the galaxies, a number of the 
galaxies were close enough to resolve a significant pop- 
ulation of stars and it is the resolved stars that we fo- 
cus on here. Previous studies of resolved stellar popu- 
lations with HST have shown that they can be used to 
answer a variety of questions. The ages and metallici- 
ties of the stars can be used to address i ssues of galaxy 
formation (e.g. iBul^ler et alJ 120 04: Br own et alJ [200.4 
iDolnhin et al.l l200ri': 'Saraicdini fc Van Duvnell2001 ^ and 
distances can be found using the tip of the red gi- 
ant branch (TRGB) meth od (e.g. iMendez et alJ 12002: 
iKarachentsev et aLll2008aj) . In this paper we present 
the galaxies' images, stellar photometry, and some ini- 
tial results, including TRGB distance measurements, the 
identification of a possible young dwarf galaxy near NGC 
4631, and the detection of LMC stars ~9° aw ay from the 
galax y center. In the second paper (Paper Il. lSeth et al.l 
l2^04^ we analyze the variation in the vertical distribu- 
tion of different stellar populations and present evidence 
that many of these galaxies host metal-poor thick disks. 

First, in §2 we describe our sample, present a detailed 
description of our image and photometry data reduction 
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process, and discuss the stellar models and photomet- 
ric transformations we use to analyze our data. In §3 
we present images and color magnitude diagrams of each 
galaxy and identify the visible stellar populations. We 
also identify a candidate dwarf galaxy lying just off the 
plane of NGC 4631, which appears to contain only young 
stars. Our method and results for TRGB distance deter- 
mination are discussed in §4 and conclusions follow in 
§5. 

2. SAMPLE DESCRIPTION &; DATA REDUCTION 

Our sample consists of 16 galaxies observed with the 
Advanced Camera for Surveys (ACS) as part of an HST 
Cycle 12 snapshot project (PI : de Jong, PI D: 9765). The 
sample was selected from t he iTullv fc Fish er (1988) and 
IKarachentsev et all l) 1993(1 catalogs to be edge-on and 
nearby. Of the 40 fields submitted for observation, 18 
fields were observed in 16 separate galaxies. 

The properties of each galaxy in the sample is given in 
Table 1. With the exception of NGC 891 & NGC 4565, 
the galaxies are all later than Sb and have maximum cir- 
cular velocities below 200 km/sec. Thus, most have a sig- 
nificantly lower mass than the Milky Way. Based on pre- 
viously published distances (m-Mut in Table 1), mostly 
determined using the TuUy-Fisher method, the nearest 
of the galaxies is NGC 55 (1.8 Mpc l|Karachentsev et alJ 
2003 aj)) and the mos t distant is NGC 4183 (19 Mpc 
ijTuIlv fc Piercell2000t l). Assuming these distances, the 
absolute B band magnitudes of the galaxies range from 
roughly -16 for NGC 4144 & NGC 5023 to roughly -20.5 
for IRAS 06070-6147 & NGC 4565. For comparison, the 
Large Magellanic Cloud has a absolute B magnitude of 
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TABLE 1 
Galaxy Sample Properties 



Galaxy 


RA(J2000) 


Dec(J2000) 


Pos. Ang. 


b [deg]'' 


Type 


Bt 


Ktot'' 


Vmax 




m-Miit 


IC 2233 


08 13 58.8 


-1-45 44 31 


-7.6 


33.06 


SBcd 


13.36 


10.75 


79 


558 


29.96' 


IC 5052 


20 52 02.9 


-69 11 45 


-38.0 


-35.80 


SBcd 


11.67 


8.88 


79 


598 


28. 80S 


IRAS 06070-6147'= 


06 07 29.9 


-61 48 26 


42.7 


-28.82 


Sc 


10.72 


8.98 


133 


1227 


31. 22'^ 


IRAS 07568-4942^ 


07 58 15.3 


-49 51 13 


-26.5 


-10.58 


SBc 


12.36 


8.04 


142 


1119 


30.068 


IRAS 09312-3248° 


09 33 21.2 


-33 02 00 


89.7 


13.64 


Sc 


12.71 


8.88 


98 


929 


30.61'^ 


NGC 55 


00 14 54.4 


-39 11 59 


-70.0 


-75.73 


SBm 


9.58 


6.25 


67 


131 


26.28' 


NGC 891 


02 22 33.0 


-1-42 20 53 


22.8 


-17.41 


Sb 


10.83 


5.94 


214 


530 


29.6lj 


NGC 3501 


11 02 47.3 


-1-17 59 22 


28.1 


63.37 


Sc 


13.67 


9.41 


136 


1132 


31. 04^ 


NGC 4144 


12 09 58.3 


-1-46 27 28 


-77.6 


69.01 


SBc 


12.16 


9.39 


67 


267 


27.93S 


NGC 4183 


12 13 16.8 


-f 43 41 53 


-14.0 


71.73 


Sc 


12.94 


9.85 


103 


931 


31.35"^ 


NGC 4206 


12 15 16.8 


-f 13 01 27 


1.0 


73.55 


Sbc 


12.83 


9.39 


124 


712 


31. 29^ 


NGC 4244 


12 17 29.5 


-1-37 48 28 


47.2 


77.16 


Sc 


11.12 


7.72 


93 


244 


28.26™ 


NGC 4517 


12 32 45.6 


-1-00 06 56 


82.6 


62.62 


Sc 


11.12 


7.33 


137 


1126 


29. 30*^ 


NGC 4565 


12 36 20.7 


-1-25 59 15 


-44.8 


86.44 


Sb 


10.58 


6.06 


227 


1228 


31.05™ 


NGC 4631 


12 42 07.7 


-1-32 32 33 


82.6 


84.22 


SBcd 


9.67 


6.47 


131 


599 


28. 18*" 


NGC 5023 


13 12 11.7 


-1-44 02 17 


27.9 


72.58 


Sc 


12.81 


8.74 


77 


406 


28. 75^ 



Note. — RA, Dec and Position Angle from K band galaxy fits described in Appendix. All other data from HYPERLEDA/LEDA iPaturel et al.l 
119951 12003) except whore noted, (a) Galactic latitude, (b) Total K-band magnitude from ljarrett et al .'f200a. (c) IRAS 06070-6147 also known as 
ESQ 1 21- G6, (d) IRA S 07568-4942 also known as ESQ 209-G9, (c) IRAS 060 70-6147 also know n as ESQ 373-G8 or UG CA 168, (f ) BottinoUi ot ^ 
1198511. (•gl lBottinel li et alj il988il. (h) assuming Ho=70 k m/s ec. fil l K arachents ev et alj I2003ar . (\) TTom:v et all I200'l^ . (kl ITullv &: Picrcg lj20odr 
fn iGavazzie^frF99Sf ). f mT FKarachentsev et all 12003bll . (n)|j ensen et al .1 12003H 



-17.6, using the Bt m agnitude from the LEDA database 
ijPatiirel et aLlllflQIiD and assuming a distance modulus 
of 18.50 l)Alve£ll2004|) . We also performed model fits 
to 2MASS K band images to obtain a consistent set of 
structural properties for use in our analysis. These fits 
are described in detail in the Appendix. In Table 2 we 
show the resulting scale length (/i/?), scale height (zq), 
and half-light height zi/2- Note that for the model used 
(Eq. 1 of the Appendix), zi/2=0.549zo. 

The HST imaging was obtained using the Wide Field 
Camera (WFC) camera of ACS. A total of 18 fields in 
16 galaxies were observed with the F606W and F814W 
filters (see Table 3). Each field was observed with one of 
two observing plans: 

• 2x2: 14 fields have two exposures of 338 seconds 
in F606W and two exposures of 350 seconds in 
F814W. 

• 1x2: 4 fields have one 400 second exposure in 
F606W and two 338 second exposures in F814W. 

In all cases where two exposures were taken, they were 
dithered to fill in the gap between the two ACS/ WFC 
CCDs. 

We developed the following data reduction pipeline for 
these galaxies: 

1. Removal of cosmic r ays using the Laplaci an edge- 
detection routine of Ivan DokkumI l)2001j) on the 
raw FLT images obtained from the HST/MAST 
archive. 

2. Sky subtraction, drizzling and further cosmic ray 
removal using the PyRAF 'multidrizzle' task (code 
written by Anton Koekemoer). 

3. Location and aperture phot ometry of obje cts in the 
image using DAOPHOT II llStetsonlll987l) . 



4. Determination of the point-spread function (PSF) 
and PSF photometry also using DAOPHOT II. 

5. Artificial star tests to obtain the completeness as a 
function of position in our images. 

All 18 fields were reduced using this same pipeline. 
Noting the lack of detailed descriptions of ACS data re- 
duction and photometry in the literature, we will now 
describe each of these steps in detail. 

2.1. Cosmic Ray Removal & Drizzling 

Based on ACS report 99-09 (Mutchler et al.), we ex- 
pect to find that approximately 1.5% of pixels (^300,000 
pixels in our final 4350x4350 images) are affected by cos- 
mic rays. Because our data set includes only two images, 
we hoped to be able to reject cosmic rays in each single 
image before combination. 

Tests of the Laplacian edge detection method on 
WFPC2 data found that it was an effecti ve method for 
remo ving CRs even on undersampled data l)van DokkumI 
120011) . We therefore adopted the same 'la_cosmic' IDL 
code for our data reduction. However, we found two 
problems when the method was applied to our data. 
First, while the code appeared to be very effective at 
removing cosmic rays, there were also a significant num- 
ber of obvious false detections in high brightness regions. 
Second, we found that removed cosmic rays still left small 
halos behind. To address the false detections we set very 
stringent limits for what was counted as a cosmic ray (sig- 
clip=6, objlim=8). To remove the residual halos, we ex- 
panded the cosmic ray removal to the four nearest neigh- 
bors of a detected cosmic ray. Although this undoubtedly 
removed some pixels unaffected by CRs it also eliminated 
all signs of the halos around removed CRs. In total, the 
altered 'la_cosmic' routine resulted in the removal and 
interpolation over ^150,000 pixels per science image in 
the FLT files. However it also left behind some obvious 
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TABLE 2 

K Band disk fit parameters 



Galaxy 


So 
[mag/D"] 


hR 

["] 


Zo 
["] 


hji/zo 


Zl/2 
["] 


hR 

[pc]=» 


Zo 


IC 2233 


18.42 


25.97 


6.87 


3.78 


3.77 


1236 


327 


IC 5052 


18.02 


53.87 


13.35 


4.04 


7.33 


1574 


390 


IRAS 06070-6147 


15.92 


16.36 


4.07 


4.02 


2.23 


1391 


346 


IRAS 07568-4942 


16.05 


27.80 


6.87 


4.05 


3.77 


1386 


342 


IRAS 09312-3248 


17.52 


38.24 


10.33 


3.70 


5.67 


2455 


663 


NGC 55 


17.38 


93.30 


42.52 


2.19 


23.34 


958 


437 


NGC 891 


15.59 


90.75 


11.86 


7.65 


6.51 


3676 


481 


NGC 3501 


16.20 


19.96 


4.15 


4.81 


2.28 


1562 


325 


NGC 4144 


18.01 


30.56 


12.68 


2.41 


6.96 


1103 


458 


NGC 4183 


17.47 


19.66 


7.21 


2.73 


3.96 


1775 


651 


NGC 4206 


17.24 


24.44 


7.45 


3.28 


4.09 


2146 


654 


NGC 4244 


17.91 


84.27 


22.14 


3.81 


12.15 


1783 


469 


NGC 4517 


17.02 


79.13 


14.64 


5.41 


8.04 


2779 


514 


NGC 4565 


16.02 


87.86 


13.67 


6.43 


7.50 


6908 


1075 


NGC 4631 


15.62 


35.49 


13.74 


2.58 


7.54 


1323 


512 


NGC 5023 


18.16 


39.81 


9.29 


4.29 


5.10 


1275 


298 



Note. — (a) Derived using distances from Table 1 or when available, Table 5 



TABLE 3 
Observed HST Fields 



Field 


Obs. Plan'' 


# Stars 


IC 2233 


2x2 


15697 


IC 5052 


2x2 


68636 


IRAS 06070-6147 


1x2 


15271 


IRAS 07568-4942 


1x2 


23582 


IRAS 09312-3248 


2x2 


21667 


NGC 55 


1x2 


281536 


NGC 55-DISK 


2x2 


253108 


NGC 891 


2x2 


33508 


NGC 3501 


2x2 


4887 


NGC 4144 


2x2 


60552 


NGC 4183 


2x2 


13103 


NGC 4206 


2x2 


8733 


NGC 4244 


2x2 


121238 


NGC 4517 


2x2 


49201 


NGC 4565 


2x2 


20192 


NGC 4631 


2x2 


104940 


NGC 4631-DISK 


2x2 


97656 


NGC 5023 


1x2 


42293 



Note. — (a) 2x2: denotes fields with two 338 second exposures 
in F606W and two 350 second exposures in F814W. 1x2: denotes 
fields with one 400 second exposure in F606W and two 338 second 
exposures in F814W. 



cosmic ray features. We note that when we had only one 
exposure in a filter, we used less stringent parameters 
for selecting cosmic rays, which may have resulted in the 
removal of real sources. However, these images were not 
used for source detections so the only affect of this would 
be occasionally spurious photometric results. 

The remaining cosmic rays were removed in the driz- 
zling process. We ran the PyRAF 'multidrizzle' task to 
create 4350 x 4350 pixel images. Due to our crowded 
fields, the default sky subtraction produced obvious mis- 
matches in sky brightness between chip quadrants, re- 
quiring us to use skytype='single', which applies a single 
sky value to each CCD image. We also found determin- 
ing the sky value using the mode (skystat='mode') rather 
than using the median (which is the default), gave bet- 



ter matched sky values. Despite these measures, some 
of our most crowded fields (e.g. NGC 4631) show a 
slight mismatch in sky brightness between the two CCDs. 
This should not affect the identification or photometry 
of stars in these images because both use locally deter- 
mined sky values. Comparing the 'median', 'crreject' and 
'minmed' modes of image combination, we found the de- 
fault 'median' mode surprisingly yielded the best results 
for cosmic ray removal. The 'lanczos3' kernel was used 
for creating the final images with a pixel scale of 0.05 
arcseconds. We found that the 'lanczos3' kernel gave a 
sharper image than the default 'square' kernel. Specifi- 
cally, within a two pixel aperture, stars in an image made 
with the 'lanczos3' kernel had '^5% more flux than in an 
image made with a 'square' kernel. For apertures larger 
than four pixels, fluxes in the two images were the same. 
All further analysis is done with the final, drizzled im- 
ages. 

2.2. Object Identification & Aperture Photometry 

All of the galaxies in the sample show a resolved stel- 
lar component. Individual stars were detected using the 
FIND routine in DAOPHOT II (Stetson 1987) with a 
threshold of 5(7 over the local background. The rou- 
tine initially returned anywhere from 15,000 sources in 
IRAS 06070-6147 to 400,000 sources in NGC 55. The 
final number of sources, shown in Table 3, count only 
sources detected in both filters. For the 1x2 fields, we 
detected objects only in the F814W filter and then used 
that source list for both filters. 

We performed aperture photometry using the PHOT 
routine with an aperture of 4 pixels and a sky annulus 
between 11 and 14 pixels. We then compiled the pho- 
tometry in both filters and required objects be detected 
independently in both filters for the 2x2 fields. This re- 
duced the number of objects by 40-75% from the original 
source lists. We note that there was a '^0.3 pixel posi- 
tional offset between the F606W and F814W drizzled 
images in all the observations. After correcting for this 
offset we matched objects within 1 pixel of eachother in 
the two filters. In the very small number of cases where 
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more than one source was found within one pixel, the 
brighter match was used. As noted above, in the 1x2 
fields, sources were found only in the F814W band - the 
positions from the F814W image were then used for both 
filters to determine the aperture magnitudes. 

We transformed all of our photometric data to the 
VEGAmag system using provided zeropoints (ACS ISR 
2004-8, De Marchi et al.) and aperture corrections (M. 
Sirianni, private communication). The zeropoints pro- 
vided are for an infinite aperture and the aperture cor- 
rections give the aperture correction from a 10 pixel to 
infinite aperture. We derived a correction from a 4 to 10 
pixel aperture using isolated bright stars from the 2x2 
fields. The resulting corrections showed very little (<0.01 
mag.) scatter and are shown in Table 4 along with the 
zeropoints and 10— oo aperture corrections. 

2.3. PSF photometry 

All of the galaxy images are very crowded, greatly re- 
ducing the accuracy of aperture photometry. To improve 
the photometry, we modeled the PSF in our drizzled 
fields. An adequate model for the ACS/WFC PSF must 
allow for spatial variations of the PSF over the wide field 
of view. Fortunately, DAOPHOT II has the capabil- 
ity to model the PSF as a combination of an analytical 
function and empirical sub-pixel corrections that depend 
quadratically on the position in the field. 

Initially we attempted to model the PSF in each field 
independently, but found that the derived PSFs were 
greatly limited by the small numbers of bright, isolated 
stars (typically 50) available in the large 4350x4350 pixel 
frames. Instead we assumed that the PSF was roughly 
constant over the course of our observations and con- 
structed a patchwork image in each filter using bright, 
isolated stars from the 2x2 fields. The resulting im- 
ages contained ^^700 stars and enabled more accurate 
determination of the spatial variation of the PSF. Out of 
the options given by DAOPHOT II, the best fitting an- 
alytical functions in both filters were Moffat functions, 
f(x) = 1/(1 -f (r/a^))''. In the F606W filter, a = 0.76 
and 13 = 2.5, while in F814W a = 0.63 and /3 = 1.5. 
The FWHM of these analytical functions is 0.86 pixels 
(0.043 arcsec) in F606W and 0.97 pixels (0.049 arcsec) in 
F814W. Based on a comparison between a variable and 
constant PSF, the spatial corrections in our PSF affect 
the photometry by at most 0.03 magnitudes. The correc- 
tion of PSF magnitudes to 10 pixel aperture photometry 
is given in Table 4. The accuracy of these corrections is 
~0.03 magnitudes. We treat this as a systematic error, 
and do not include it in individual errors of the stars. 
The computed PSF has a radius of 7 pixels (0.35 arcsec). 

Figure la shows the resulting PSF magnitude color- 
magnitude diagram (CMD) for the NGC 55-DISK field. 
The NGC 55-DISK field CMD is the cleanest of our 
CMDs due to the galaxy's proximity and the field's rel- 
ative lack of dust. The clear features seen in this CMD, 
which will be discussed in §3, demonstrate the quality of 
our photometry. 

2.4. Artificial Star Tests 

To assess the completeness limit as a function of posi- 
tion, magnitude, and color, we inserted large numbers of 
artificial stars in our images and then attempted to re- 
cover them with our data reduction pipeline. These stars 



were input in grids separated by 20 pixels with an overall 
random position offset. In this way, ~38,000 stars were 
added to a field, each with a random F606W VEGAmag 
between 18 and 29 and a F606W-F814W color between 
-1.0 and 3.0. We then ran the stars through a pipeline 
identical to the one outlined above for finding stars and 
determining their photometry. This was done with dif- 
ferent random position offsets to generate as many as 
four million artificial stars distributed throughout each 
field. The large computing time required to run these 
tests was enabled through use of the distributed com- 
puting system. Condor^. 

Figure 2 shows the completeness as a function of scale 
height in NGC 4244 for stars at a variety of magnitudes. 
Away from the plane of the galaxy, in a non-crowded 
area, the 50% completeness level of our data is ~27 
mag in F606W and ~26 mag in F814W. However, the 
crowding in the midplane of the galaxy drastically re- 
duces the completeness, by as much as 35% at moder- 
ate magnitudes. Note that Figure 2 contains the results 
from artificial stars of all colors, but for a specific range 
in color the completness can be higher or lower. For 
instance, the F814W 50% completeness in non-crowded 
areas at F606W-F814W between -1 and 1 is ~26.5, but 
rises steeply at redder colors. Figure 3 shows a typical 
20% completeness line as a function of color in each of 
the color-magnitude diagrams. The scatter from galaxy 
to galaxy in the completeness values is a few tenths of a 
magnitude. 

2.5. Stellar Models & Photometric Transformations 

Critical to the interpretation of our data is the abil- 
ity to connect the ACS VEGA magnitude system for the 
F606W and F814W filters with existing stellar models 
and previous data. We rely on two sources for perform- 
ing this photometric transformation. First, we used a 
set of the Padova isochrones and single stellar popula- 
tion models transformed into the ACS VEGAmag system 
using synthetic spectra and ACS/WFC instrument re- 
sponse (L. Girardi, private communicati on) . These mod- 
els are the standard lGirardi et alJ l)2000|) isochrones with 
tracks at Z=0.0001, 0.0004, 0.001, 0.004, 0.008, & 0.019 
([Fe/H]=-2.3, -1.7, -1.3, -0.7, -0.4, 0.0) with ages between 
4 Myr and 18 Gyr. 

Second, in very limited cases, we used transformations 
from the ACS VEGAmag system to Johnson V & I pro- 
vided by M. Sirianni (private communication). However, 
these transformations are potentially inaccurate as they 
attempt to infer the spectrum of a star from limited color 
information and are valid only for -0.5<V-I<2.3. 

The VEGAmag F606W-F814W color is similar to 
Johnson V-I, but has a somewhat smaller wavelength 
span, meaning that the differences in color are somewhat 
more moderate. For instance. Red Giant Branch stars 
are -0.3 magnitudes bluer in F606W-F814W [VEGA- 
mag] than in V-I. Unless explicitly stated otherwise, all 
magnitudes in the paper are in the VEGAmag system. 

3. COLOR-MAGNITUDE DIAGRAMS 

^ The Condor Software Program (Condor) was developed by 
the Condor Team at the Computer Sciences Department of the 
University of Wisconsin-Madison. All rights, title, and interest in 
Condor are owned by the Condor Team. 



A Study of Edge-On Galaxies with HST/ACS. I. 



5 



TABLE 4 

Photometric Constants 



Filter 


Zeropoint 


10— CO Ap. Cor. 


4-10 Ap. Cor. 


PSF-10 Cor. 




[mag] 


[mag] 


[mag] 


[mag] 


F606W 


26.398 


0.091±0.001 


0.079±0.006 


0.017±0.027 


F814W 


25.501 


0.090±0.001 


0.095±0.009 


0.019±0.025 




F606W-F814W [VEGAmag] F606W-FS14W [VEOAmag] 



Fig. 1. — (a) Color-Magnitude Diagram (CMD) of the 253,108 stars in the NGC 55-DISK field. Contours are used where the density of 
points becomes high. The contours are drawn at densities of 75, 100, 150, 200, 250, 350, 500, 750, 1000, and 1500 stars per 0.1 magnitude 
and color bin. (b) Cartoon CMD showing the location of Main Sequence (MS), Red Giant Branch (RGB), Asymptotic Giant Branch 
(AGB), Red Helium Burning Sequence (RHeB), and Blue Helium Burning Sequence (BHeB) each of which is described in §3. 




0.0 f .... I .1 ....... I. I 1 

-100 -50 SO 100 150 
Arcseconds from midplane 

Fig. 2. — The fraction of artificial stars detected as a function 
height above the midplane of the galaxy in NGC 4244. The lines 
give the completeness in 0.5 mag wide bins around a magnitude 
shown on the right-hand side. Dashed lines show ±3zi/2i the min- 
imum height used for the TRGB determination in §4. 



In the sample as a whole we have determined photom- 
etry for over 1.2 million objects (see Table 3), a vast 
majority of which are stars in the sample galaxies. 

Figure 3 shows images and color-magnitude diagrams 
(ClVIDs) of each of the galaxies. The CMDs show aU 



the stars detected in each field. Of all the ClVIDs in our 
sample, the NGC 55-DISK CMD reaches the faintest ab- 
solute magnitudes. In Figure lb we illustrate the distinct 
features visible in its CMD. Many of the other CMDs in 
Figure 3 share these features at fainter apparent mag- 
nitudes. These visible features are due to the following 
stellar populations: 

• Main Sequence (MS) - At the bright absolute mag- 
nitudes at which we are observing (Mf814 < 0), 
MS stars are younger than a few 100 Myr. At the 
bright end, the MS and Blue Helium Burning se- 
quence overlap. 

• Red and Blue Helium Burning Sequences (RHeB 
and BHeB) - These sequences contain young mas- 
sive stars undergoing core helium burning. Along 
both these sequences, the ages and masses of 
the stars are proportional to the luminosity, with 
younger, more massive stars populating the bright 
end. The faintest stars along the RHeB have 
an age of ~200 Myr, while the BHeB branch 
extends down to ^^600 Myr jP ohm-Pa lmcr et alJ 
I2OOI FDohm-Palmer fc Skillm£di2002jl . 

• Asymptotic Giant Branch (AGB) - The AGB con- 
tains stars of intermediate age, after completion of 
their core-helium burning stage. The observed por- 
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-10 12 3 

F606W-F81 4W [VEGAjuag] 

Fig. 3. — (left) - Combined F606W and F814W images of each field. Dashed fines indicate 3z]^/2 cross indicates galaxy center as 
determined from fits to K band 2MASS data (see Appendix). The images are 218"x218". The box shown in the NGC 4631-DISK field 
gives the location of image shown in Fig. 4. Note that color images of these pictures are presented in the electronic version, (right) — Color 
magnitude diagrams of all the objects in each field. Contours are shown where densities exceed 75 stars per 0.1 magnitude and color bin. 
Contours have the same levels as those shown in Fig. 1. The dashed line indicates the typical 20% completeness limit. The solid horizontal 
line (when present) indicates the TRGB. The nearly vertical feature at F606W-F814W~0.4 in IRAS 06070-6147 is due to LMC field stars. 
Distance moduli m-Mnt or m-MTRGB come from Tables 1 and 5. 
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3. — continued 



tion of the AGB is brighter and redder than the 
RGB. Based on the Girardi isochrones, the AGB 
stars observed here mostly have metahicities be- 
tween [Fe/H]=-0.7 and 0.0 and ages between 100 
Myr and a few Gyr. At absolute magnitudes fainter 
than the red giant branch tip Mf814w ^ —4 the 
AGB stars cannot be distinguished from the RGB 
stars 



• Red Giant Branch (RGB) - The RGB is domi- 
nated by old stars undergoing shell-hydrogen burn- 
ing. These stars have ages > 1 Gyr and represent 
the oldest of the components visible in our CMDs. 
The exact position of the RGB depends primar- 
ily on metallicity, with more metal-rich stars being 
redder and fainter in our CMDs. 
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Fig. 3. — continued 



Inspection of Figure 3 shows that most of these com- 
ponents are clearly visible in IC 2233, IC 5052, NGC 55 
(both fields), NGC 4144, NGC 4244, NGC 4631-DISK, 
and NGC 5023. These CMDs clearly show evidence for 
both young and old stellar populations. The effects of 
dust are also obvious when comparing the color mag- 
nitude diagrams of the central NGC 55 and NGC 4631 
fields with those taken further out in the disk (NGC 55- 



DISK and NGC 4631-DISK). In both cases the central 
CMDs show significant redward blurring of the MS and 
RHeB, both of which are more clearly detected in the 
less obscured '-DISK' fields. 

We note that a vast majority of stars detected in each 
field appear to be in the target galaxies. This can be seen 
in the clarity of their CMDs and from the spatial distri- 
bution of the stars (see Paper II). A significant popula- 
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Fig. 3. — continued 



tion of Galactic foreground stars are visible in the CMDs 
of the three galaxies with the lowest galactic latitudes 
(see Table 1), IRAS 07568-4942, IRAS 09312-3248 and 
NGC 891. They are seen most prominently at F606W- 
F814W > 0.5 and F814W < 23. 

Before discussing distance measurements using the 
RGB tip in §4, we briefly discuss two unusual features 
seen in the CMDs of IRAS 06070-6147 and NGC 4631- 



DISK. All further interpretation of the CMDs wiU be 
presented in Paper II, in which we examine the verti- 
cal distribution of the different stellar populations and 
find evidence for old and metal poor thick-disks in these 
galaxies. 

3.1. A Detection of the Large Magellanic Cloud (LMC) 
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Fig. 3. — continued 



One unusual feature is found in the CMD of IRAS 
06070-6147. In the CMD, a clear locus of points can be 
seen extending between 20th and 24th magnitude and 
at a color of ~0.5. These are field LMC stars found ^9 
degrees (8 kpc) from the LMC cent er, near the outermost 
clusters associated with the LMC l)Irwinlll991|) . 

These sta rs fit well to isochrones with a distance mod- 
ulus of 18.5 llAlvesl2004D . an age of a few Gyr, and metal- 



licities of [Fe/II]=-0.4 to 0.0. This provides a serendip- 
tous confirmation of the accuracy of our photometric cal- 
ibration and isochrone band pass transformation. It also 
agrees with recent work bv IGallart et al.l lj2004j) which 
suggests that the extended 'halo' of the LMC is made up 
of stars as young as ~2.5 Gyr. 

3.2. A Young Dwarf Galaxy Candidate in NGC 4631 
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Fig. 4. — An F606W image showing a close-up of the section of 
NGC 4631-DISK shown as a box in Figure 3. The circle has a radius 
of 4.5"and is centered at o = 12h41m50.30s and <5 =32°31'1.9". 
The X marks the location of previously identified galaxy NGP9 
F268-1993301 lOdewahn fc AlderineHQijai . A clear overdensity of 
stars is seen inside the circle. 
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Fig. 5. — (left) - A CMD of the stars inside the circle shown in 
Figure 4. Overplotted is an isochrone with an age of 30 Myr and 
[Fe/H]=-0.7. (right) - The CMD of a comparison field located at 
the same galactic latitude as the dwarf, but offset a short distance 
in the longitudinal direction. An absolute magnitude scale is shown 
assuming m-M=29.42. 



Inspection of the NGC 4631-DISK field image showed 
a small object marked by a grouping of blue stars sitting 
significantly off the plane of the disk. This object is 
shown in Figure 4 and is centered at a =12h41m50.30s 
and 6 =32°31'1.9". The object seems to be contained 
within a radius of 4.5"as shown in Fig. 4. Assuming 
a distance modulus of 29.42 to the object, this gives a 
radius of '^170 pc. The object is located ~2 kpc above 
the plane of NGC 4631. Based on its size, which is too 
large to be a star cluster, we identify this object as a 
dwarf galaxy candidate, NGC 4631 DwA (dwarf A). 

This dwarf may have be en previously discovered by 
lOdewahn fc Alderind l)1995D with the name NGP9 F268- 
1993301 using POSS-I plates. However, their position 
(a=12h41m49.89s (5-|-32°31'00.5", shown with an X in 
Fig. 4) and red color (0-E~2) suggests that they suffered 
from confusion with the bright red star (most likely a 



foreground object) seen to the lower left of the circle in 
Figure 4. 

The left panel of Figure 5 shows the color-magnitude 
diagram of the stars in the circle shown in Figure 4. For 
comparison, we show a field CMD at the same height 
above the plane of NGC 4631 but offset in the longitudi- 
nal direction. NGC 4631 DwA shows a clear overdensity 
of blue stars. To show this quantitatively, we compare 
the stars in the dwarf to those within the same range 
of height above the plane, but lying outside the dwarf, 
yielding an area '^18 times the area of the dwarf. We 
find 26 stars bluer than F606W-F814W=0.3 associated 
with the dwarf, but only 15 in the entire area outside 
the dwarf. This gives an expected rate ~0.8 stars bluer 
than F606W-F814W=0.3 in the field of the dwarf, imply- 
ing that the dwarf has a >30 times overdensity of blue 
stars. By contrast, the expected rate for stars redder 
than F606W-F814W of 0.3 in the field of the dwarf is 
~50 stars, while the dwarf contains only 47. This sug- 
gests that if any old population exists, it is small enough 
to contain few stars near the tip of the RGB. The good fit 
of the isochrone shown in Figure 5 to the blue stars and 
the bright red stars (probably RHeB stars) in the dwarf 
suggests an age of ~30 Myr. We found that isochrones 
with a wide range of metallicities and ages between 10- 
50 Myr can also fit the data quite well. These tests also 
suggest that the dwarf is located at the same distance 
modulus (assumed to be 29.42) as NGC 4631. 

We can also calculate the integrated magnitudes of 
NGC 4631 DwA to compare to single stellar population 
models. The total fiux from the dwarf is measured by 
aperture photometry using an aperture of 4.5" and a sky 
annulus between 5" and 6". The contribution from NGC 
4631 's field stars is then determined by finding the fiux 
using the same aperture and sky annulus at 100 ran- 
dom positions with the same height above the plane, 
but avoiding NGC 4631 DwA and the bright star vis- 
ible just to it's lower left in Figure 4. The mean of 
these we use as a background and the standard devi- 
ation is the error in this background. The final fiux 
for the dwarf is determined by subtracting off the back- 
ground and the final error combines the dwarf's error 
in the flux with the error in the background. Using 
this error and the background subtracted flux we find 
that the F606W magnitude is 20.7±0.2 (a 5a detection) 
while the F814W magnitude is 20.3±0.4 (a 3cr detection). 
These magnitudes correspond to surface brightnesses of 
25.2 mag arcsec"^ in F606W and 24.8 mag arcsec"^ in 
F814W. Assuming an age of 30 Myr and [Fe/II]=-0.7, 
we find a mass of 2.2±0.4xl0^ Mq from single the stel- 
lar population models (see §2.5). We note that despite 
being rather uncertain, the F606W-F814W color of 0.4 
matches the expected value from this single stellar popu- 
lation, suggesting that any existing older population does 
not make a significant contribution to the total luminos- 
ity of NGC 4631 DwA. The dwarf's absolute magnitude, 
Mf606W = —8.7 is only slightly brighter than t he recently 
disco vered 'faintest galaxy' Andromeda IX (Zucker et alJ 
I200I . 

The properties of NGC 4631 DwA - a radius of 170 
pc, an age of ~30 Myr, and a mass of ~2xl0'* M0 are 
most similar to an OB association. However its location 
~2 kpc above the disk makes it unusual. If the object 
is indeed bound it would have a circular velocity of ~7 
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km/sec. However given that the only stars we see are 
young, it is possible that it is an unbound and gradually 
dissolving object. 

One possible explanation for the creation of this ob- 
ject is that it is caused by a condensation of gas inter- 
acting with NGC 4631. Figures 4 and 5 of Rand (1994) 
show five HI spurs surrounding NGC 4631. Although 
none of their positions overlap with the position of the 
dwarf, two spurs meet up with the disk at a position 
just opposite the dwarf. This explanation for the for- 
mation of the object would make it akin to young tidal 
dwarf galaxies and clu sters observed in the tails of in- 
teracting gala xies (e.g. lSchweizeI^ 119781 iDuc & Mirabell 
^98; Tempori n et al.l l2'b03). Higher resolution HI obser- 
vations would make it possible to determine if the object 
is affiliated with a tidal tail. 

4. TRGB DISTANCES 

The tip of the red giant branch (TRGB) method can 
be used to determine distances to objects with an old 
and metal poor ([Fe/H] < -0.7) resolved stellar popu- 
lation. Of the sample of galaxies presented here, eight 
fields in six galaxies were close enough to clearly distin- 
guish an RGB tip. To determine the TRGB magnitude 
we use the discontinuity detection method as presented 
in Mendez et al. ( 2002). 

4.1. Star selection and method 

Before considering the method by which we determine 
the TRGB, we need to refine the sample of stars. To pre- 
vent our result from being affected by in situ dust and 
young stellar populations and to avoid incompleteness 
(Fig. 2), we consider stars only above 3|zi/2|. This height 
is above >90% of the K band light on the 2MASS images 
and is above any obvious dust features on our images, 
except in NGC 4631 which is known to have a lar ge ex- 
traplanar dust component l|Howk fc Savagelll999D . The 
CMDs constructed from these stars show the presence of 
very few young stars. This can be seen in Figure 6, which 
shows the CMD of stars above 3|zi/2l except in NGC 
4631 and NGC 4631-DISK fields where we used a lower 
limit of 7|z]^/2l in order to better avoid the dust, crowding 
and young stellar populations seen at lower scale heights. 
For the NGC 4631 field, we also used an upper limit of 
10|z]^/2l to avoid stars in the nearby companion, NGC 
4627. 

Comparison of the CMDs in Fig. 6 to the full field 
CMD in Figure 3 shows the lack of young RHeB stars and 
the greater clarity of the RGB due to a lack of crowding. 
After selecting these stars, we correct their magnitudes 
for foreground extinction using Schlegcl et al. ( 1998) val- 
ues for E(B-V) and assuming Af814w=2.716xE(B-V) 
and Af606W=1-796xE(B-V). Finally, we make a color 
cut (0.3 < F606W-F814W < 1.6) to eliminate contami- 
nation of the upper main sequence and redder AGB stars. 

Figure 7 graphically shows the method used to de- 
termine TRGB magnitude and error in each field. Us- 
ing the magnitudes and errors of the selected stars, we 
construct a Gaussian-smoothed luminosity function, to 
which each star contibutes a unit area Gaussian with 
th e width determin ed by the error (see Appendix A 
of ISakai et al.l ir996(l . This is shown in the top pan- 
els of Figure 7. The tip of the red giant branch can 
be seen as a sudden increase in the luminosity function 



at a given F814W magnitude. This magnitude can be 
quantitatively determined using an edge-detection algo- 
rithm. We use the logarithmic edge-detection algorithm 
presented in^endcz et al. (2002). The response of this 
edge-detection algorithm is shown in the middle panels 
of Figure 7. Finally, random errors in the TRGB magni- 
tude are determined using a Monte Carlo style simulation 
where the TRGB magnitude is determined 500 times for 
stellar magnitudes which have been randomly resamplcd 
with replacement. We also include the magnitude errors 
by adding a Gaussian random error scaled by each stars 
photometric error. The width of a Gaussian fit to the 
determined TRGB magnitudes gives us th e uncertainty 
of the derived TRGB value (Me ndez et 311120021) . A his- 
togram of these Monte Carlo tests can be seen in the 
bottom panels of Figure 7. 

4.2. TRGB magnitudes 

The results of the TRGB determination for the eight 
fields in which a clear RGB could be distinguished are 
shown in Table 5 and in the CMDs in Fig. 6. The F814W 
TRGB magnitudes range from 22.54 for NGC 55 to 25.34 
for NGC 4631. The uncertainties are typically ~0.03 
magnitudes, comparable to the err ors found using a simi- 
lar method in'Men dez et alJ l)2002t ) and to the systematic 
errors in our magnitudes (§2.3). In all cases, the value 
derived for the TRGB fell within one sigma of the Gaus- 
sian fitted to the Monte Carlo tests. The second column 
of Table 5 shows the number of stars used in the TRGB 
determination of each galaxy. 

For each of the galaxies shown in Table 5, a clear 
TRGB is seen (Fig. 6). For the other galaxies, we are 
limited in our ability to determine TRGB distances by 
incompletness. The F814W 90% completeness level is 
^25.5 in uncrowded areas at RGB colors and drops off 
quickly towards fainter magnitudes. The only two other 
possible detections are NGC 891 and NGC 4517, both of 
which we measured to have TRGB magnitudes of 25.60. 
However, examination of the CMDs reveals that it is im- 
possible to exclude the possibility that the stars below 
that magnitude are AGB or RHeB stars. In addition, 
the Monte Carlo tests find a wide range of values for 
the TRGB and their results are not consistent with the 
a TRGB magnitude of 25.60. We suggest that the de- 
rived TRGB magnitude for these galaxies is a lower limit. 
For the other 8 fields in our sample, there were typically 
too few stars above 3|zi/2| to distinguish any red giant 
branch. 

4.3. Distance Moduli 

The I-band absolute magnitude of the TRGB has 
been found to be at Mj ^ —4 via many methods 
(e.g.lMendez et al.ll200l iLee et al.lll993HFerrarese et alJ 
12000 ). Because the photometric system we use is similar 
to the Johnson system, we can expect that the TRGB is 
at a similar F814W absolute magnitude (Mf814w)- In or- 
der to accurately determine the distance, we need to (a) 
determine the observed TRGB magnitude in the John- 
son I system, or (b) determine the exact value Mf814W 
of the TRGB. We take both of these approaches below. 

Using the transformations provided to us by Siri- 
anni, we transformed the stellar data for each galaxy 
to the V and I filter and redetermined the TRGB ex- 
actly the same as above, but using V-I cuts of 0.5 and 
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Fig. 6. — The CMDs of stars used in TRGB determination. The hne in each pane indicates the magnitude found for the TRGB. Stars 
were selected to have 0.3<F606W-F814W<1.6, and to be well above the plane of the galaxy (>3zx/2 fo'' s-U galaxies other than NGC 4631). 
These CMDs show a much greater clarity than those in Fig. 3 because of the reduced effects of dust and crowding. 



TABLE 5 

TRGB RESULTS 



Field 


# of stars used 


F814Wtrgb 
[VEGAmag] 


Itrgb 
[Johnson] 


(F606W-F814W)_3 5 
[VEGAmag] 


(V-I)_3.5 
[Johnson] 


[Fe/H] 


m-MTRGB 


IC 5052 


18974 


24.82±0.03 


24.84±0.04 


1.15 


1.48 


-1.22 


28.90 


NGC 55 


53984 


22.55±0.03 


22.55±0.03 


1.08 


1.42 


-1.41 


26.63 


NGC 55-DISK 


56651 


22.54±0.03 


22.55±0.03 


1.04 


1.36 


-1.62 


26.62 


NGC 4144 


19132 


25.28±0.02 


25.29±0.03 


1.06 


1.39 


-1.50 


29.36 


NGC 4244 


27797 


24.12±0.03 


24.14±0.03 


1.08 


1.41 


-1.45 


28.20 


NGC 4631 


11204 


25.34±0.03 


25.37±0.03 


1.13 


1.47 


-1.25 


29.42 


NGC 4631-DISK 


12058 


25.33±0.03 


25.36±0.04 


1.06 


1.38 


-1.54 


29.41 


NGC 5023 


8376 


25.02±0.03 


25.04±0.03 


1.03 


1.33 


-1.71 


29.10 



1.9. The results are shown in Table 5. By calibrat- 
ing the TRGB with Cephcid distances, iFerrarese et al.l 
(pOOO), found the TRGB to be at IVli = -4.06 ± 0.07 
(random) ±0.13 (systematic). The magnitudes of Itrgb 



agree with F814Wtrgb to within 0.03 magnitudes for 
each of the 8 fields, but on average Itrgb is fainter by 
^0.02 magnitudes. This suggests that the TRGB is at 
Mf814W 4.08. 



A Study of Edge-On Galaxies with HST/ACS. I. 



15 





20.5 21,0 21,5 2Z0 22,5 23.0 23.5 24.0 21 22 23 24 25 26 

F606W-F814W [VEGAraag] F606W-F814W [VEGAmag] 



Fig. 7. — TRGB determination for eight fields in six galaxies. In each plot, the top panel shows the Gaussian-smoothed luminosity 
function as a function of the F814W magnitude, the middle panel shows the edge detection response to the luminosity function and the 
bottom panel shows a histogram of the edge-detection peak for 500 monte carlo simulations. The dashed line gives the position of the 
determined TRGB value. For more details see §4.1. 



Alternately, we can use the theoretical Padova 
isochrones (see §2.5), which exist in both filter systems. 
For both sets of isochrones we take all tracks older than 
2 Gyr and with [Fe/H] between -2.3 and -0.7. Figure 8 
plots the cumulative distribution of the TRGB magni- 
tude in both the F814W VEGAmag system (sohd) and in 
Johnson I (dashed). A great majority of the tracks have 
TRGB values in both systems between -3.8 and -4.1. The 
mean, median and brightest values of the absolute mag- 
nitude of the TRGB are all ^^0.03 magnitudes brighter 
in F814W VEGAmag than in Johnson I, in good agree- 
ment with the empirical calibration above. We therefore 
assume that Mf814w = —4.08 for determining distances 
from the observed TRGB magnitudes. The resulting dis- 
tance moduli are shown in the final column of Table 5. 
The errors in the distances are the random error asso- 
ciated with the TRGB measurement a dded in quadra- 
ture with the 0.07 mag error from the iFerrarese et al.l 



l)200nt) calibration, resulting in a ^^0.08 mag random er- 
ror. The systematic error is ap proximately 0.1 5 mag- 
nitudes, slightly higher than the IFerrarese et aD l)2000f) 
value due to the change in filter system. 

Comparing these distances to the literature values 
compiled in Table 1 reveals good agreement on two galax- 
ies, with less good agreement for the other four. We note 
that our limit at F814Wtrgb of ~25.5 gives a distance 
modulus limit of ~29.6. Based on previous distance 
measurements, we do not expect to detect the TRGB 
in any other galaxies other than those listed in Table 5 
except possibly for NGC 891 and NGC 4517, whose non- 
detection we discussed above. 

We now compare our distances with previous distances 
on a case-by-case basis. The listed A(m-M) value is 
found by taking our distance modulus and subtracting 
the previous measurement. 
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Fig. 8. — Cumulative distribution of Padova isochrones with 
TRGB magnitudes brighter than a given absolute magnitude. Solid 
line shows the distribution for MFgi4w VEGAmag isochrones, 
while the dashed line shows the Johnson Mi values. 



• IC 5052 (A(ni-M)=0.10): The previous distance 
measurement is from the IR TuUy-F isher rela- 
tion as given in jB ottincUi et alJ l)1988D . In fact, 
iBottinelli et aLl lfl988ll gives two values for the dis- 
tance modulus, calculated using different TuUy- 
Fisher relations, the mean of which gives the dis- 
tance modulus we have found. 

• NGC 55 (A(ni- M)=0.35): The distance gi ven in 
Table 1 is from iKarachentsev et alJ l)2003a|) , who 

reinterprets older HI TuUy-Fisher data collected 
from various sources by Puchc & Carignan (198^ 
to determine the distance of 1.8 Mpc (m-M=26.28). 
Previous measurements using t he magnitudes of su- 
pergiants ijPritchet et aUlToS^ from ground-based 
narrow band imaging find an even closer distance of 
1.38 Mpc (m-M=25.7) which disagrees even more 
severely with our measured distance. We suggest 
that these previous measurements may be unreli- 
able, especially the TuUy-Fisher relationship in this 
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irregular galaxy. We also note that our distance 
modulus agrees within the erro rs to the TRGB & : 
Cepheid distance to NGG 300 ijButler et al.ll2004l: 
[Precdman et al. 2 001|) . another galaxy in the Sculp- 
tor Group. Assuming similar distance moduli, 
the separation of 8 degrees between these galaxies 
translates to a physical distance of '^300 kpc. 

• NGC 4144 (A(m-M)^1.4 3): The Table 1 m- 
M=27.93 value is from the iBottinelH et all l|1988l^ 
IR TuUy-Fisher relation. Other estimates give 
wildly varying distan ces: a previous B band 
TuUy-Fisher distance ijBottinelh et alJlT985^ gave 
m- M=28.87, while a brightest s t ar me asurement 
by iKarachentsev fc Drozdovskvl l|1998fl gives m- 
M=29.9. Our distance measurement falls between 
these latter two. 

• NGC 4244 (A(m-M)=-0.06): The close agreement 
between our distance and the distance in Table 1 is 
unsurprising because both were derived using the 
TRGB method from HST data. We note that a 
previous B-band TuUyj-Fishcr measurement finds 
m-M=27.77 (|Bot'tinelii et al]ll984il . ~0.5 magni- 
tudes away from the TRGB distance. 

• NGC 4631 (A(m-M)=1.24): The previous mea- 
surement of m-M=28.1 8 is from the B-ban d 
TuUy-Fisher relation of iBottineUi et all ljl985j) . 
Measurements of the H I and CO TuUy-Fisher 
ijSchoniger & Sofu'3ll994D relation suggest an even 
closer distance modulus of 27.5 to 27.8, in even 
worse agreement with the TRGB distance. We note 
that this galaxy is undergoing interactions l|R.andl 
^994) and has nearby companions, both of which 
could make the TuUy-Fisher distances less reliable. 

• NGC 5023 (A(m-M)=0.35): This difference in m- 
M is within the very large error (0.5 mag) given by 
[BottineUi et al. (1985) for the B-band TuUy-Fisher 
measurement of this galaxy. 

To provide another check on our distance measure- 
ments, we compare our TRGB distances to the K 
band/Vmax TuUy-Fisher relation presented in IVerheiienI 
lf2f)(m : 

Mk = -1.62 ± 0.49 - 8.5 ± 0.2 log(2Vmax) (1) 

Figure 9 shows the comparison between this TuUy-Fisher 
relation and the distances we have derived. The magni- 
tudes shown were calculated using the 'total' Kg mag- 
nitude (Ktot in Table 1) from Jarrett et al. ( 2003) and 
our TRGB distances. The Vma x values are those shown 
in Table 1 ijPaturel et al.ll2003|) and represent the max- 
imum circular velocity in that galaxy. We note that K' 
filter used for deriving the TuUy-Fisher relationship has 
a similar bandpass to the Kg filter used in 2MASS, and 
so any differences between the magnitude systems should 
be small. For the six galaxies that we have derived dis- 
tances to, the scatter around the TuUy-Fisher relation is 
0.65 magnitudes. This is similar to the 0.71 magnitude 
scatter around the relation found by Verheiicn (2001). 

Overall, we find that there are plausible reasons for 
the differences between our distances and those previ- 
ously measured. The agreement found between different 
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Fig. 9. — A comparison of the derived TRGB distances (data 
points) vs. the K band TuUy-Fisher relation from Verhciien (200^ 
(Une). In addition to the six galaxies with detected distances, we 
show NGC 891 and NGC 4517 assuming a lower limit to the dis- 
tance modulus of 29.68. 

fields of the same galaxy as well as the consistency of 
our TRGB distances with the K band TuUy-Fisher re- 
lationship and the previous TRGB measurement (NGC 
4244) all give us reason to believe the distances given 
here are more accurate than previous distances found in 
the literature. We also note that the relative distances 
between the galaxies should be accurate to within the 
~0.03 magnitude errors in the TRGB determination. 

4.4. Metallicity of the RGB 

Table 5 shows the color of the RGB at an absolute 
magnitude of -3.5 in both the Johnson system, (V-I)_3 5 
and the ACS VEGAmag system, (F606W-F814W)_3.5- 
These values were found by fitting a Gaussian to the dis- 
tribution of stars within an 0.2 magnitude band centered 
at an absolute magnitude of -3.5 (assuming TRGB at 
Mf814W=-4.08 an d Mt=-4.06) . T he value for (V-I)„3.5 
has been used bv lLee et al.l () 1993^1 to derive the metal- 
licity as follows: 

[Fe/H] = -12.64 + 12.6(V - I)_3.5 - 3.3(V - 5 

(2) 

The [Fe/H] values from this equation are shown in Ta- 
ble 5 and range from -1.2 to -1.7. Comparison of 
the (F606W-F814W)_3.5 values (Table 1) with Padova 
isochrones shows similar results, with the data falling 
just to the metal rich side of the [Fe/H]=-1.3 isochrone. 
This confirms that the population of stars we are using 
for TRGB distance determination are metal poor and 
suggests that they are part of some kind of halo or thick 
disk system. Note also that the outer fields of NGC 55 
and NGC 4631 have lower metallicities than the inner 
fields as expected in some disk formation models. We 
will explore these trends more fully in Paper II. 

5. SUMMARY 

In this paper we have presented the basic properties 
and data reduction of a sample of edge-on galaxies ob- 
served with ACS camera on HST. We summarize our 
findings below: 

• A number of the sample galaxies are close enough 
to be partially resolved into stars with an identifi- 
able Red Giant Branch, Asymptotic Giant Branch, 
Main Sequence, and Helium Burning sequences. 
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• We identify a candidate dwarf galaxy near NGC 
4631 that has a mass, age, and size similar to an 
OB association, but lies 2 kpc above the disk. We 
also find LMC stars in the IRAS 06070-6147 field, 
~9° from the LMC center. 

• We determine scale lengths and heights for each 
galaxy using single-disk fits to K-band 2MASS im- 
ages. 

• TRGB distances are derived to eight fields in six 
galaxies. These are found to have distance moduli 
between 26.62 and 29.42. 

• The metallicities of the high latitude stars appear 
to be metal-poor with [Fe/H] < -1. 

This work is the first in a series of papers, the next of 
which focuses on the stellar populations of the galaxies 
as a function of scale height (Paper II). Future work in- 



cludes a third paper on the star cluster content of the 
sample galaxies. 

The authors would like to acknowledge the help of Tom 
Jarrett in obtaining 2MASS data, Thomas Brown and 
Marco Sirianni for their help in obtaining ACS photomet- 
ric calibrations and Leo Girardi for providing isochrones 
in the ACS filter system. Additionally we'd like to thank 
Simone Bianchi, Paul Hodge, Kevin Covey, Vandana De- 
sai, Peter Yoachim, and our referee for helpful discussions 
and comments. Julianne J. Dalcanton was partially sup- 
ported by the Alfred P. Sloan research foundation. This 
research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space 
Administration. This work was supported by HST-GO- 
09765. 



APPENDIX 
K BAND DISK FITS 

We describe here the fits performed on 2MASS K band images to determine a consistent set of structural properties 
in our galaxies. We obtained 2MASS K band images of each galaxy in our sample from the Large Galaxy Atlas 
(parrctt ct al. 2003). Only one of the galaxies, NGC 3501, was not in the Large Galaxy Atlas, and this we obtained 
via the 2MASS Extended Source Image server^. 

We fit these K band images to an edge-on single disk model from Ivan der Kruit & Sear lei l|1981fl : 

I](i?, z) = T.o{R/h)Ki{R/hR)sech'^{z/ zo) (Al) 

where R and z are the radial and vertical coordinates. So is the edge-on central surface brightness, h]^ is the radial 
scale length and zq is the vertical scale height of an isothermal stellar population. Note that zq is different from an 
exponential scale height. For this reason we use the zi/2 parameter, which gives the height above and below the plane 
that contains 50% of the K band light to facilitate comparison with other studies. For our adopted profile in Eq. 1, 
Zi/2=0.549zo. 

Bef ore fitting, the 2MASS images were cropped and sources were identified using SExtractor ijBertin fc Arnoutl 
11996ft . Sources other than the target galaxy were masked and the masked image was fit to the single disk model of 
Eq. 1 using a Levenberg-Marquardt least-squares fit with a uniform weighting of all the unmasked pixels. For galaxies 
with obvious bulge components (IRAS 06070-6147, NGC 891, NGC 3501, NGC 4206, NGC 4244, and NGC 4517), 
we masked out the central portion of the galaxy with a radius between 5 and 50 arcsec depending on the size of the 
bulge. We attempted to fit a second disk component to the galaxies, but found that it did not increase the quality 
of the fit. Four galaxies, NGC 55, NGC 4183 and NGC 4631, showed significant asymmetry complicating the fitting 
process, while one galaxy, NGC 891 had a significant dust lane even at K band. The formal errors on zq and hji were 
at most 4%, however, the actual errors will in most cases be dominated by source structure in the galaxies. 

The results of the single disk fits are shown in Table 1 (RA, Dec and position angle, given E of N) and Table 2 (Sq, 
hfi, and zq). Also shown in Table 2 are the ratio of scale length to height /i^j/zo, the half-light height zi/2 and the hn 
and Zq values in parsecs assuming the distances in Tables 1 & 5 (where available, TRGB distances from Table 5 are 
used). The values for the scale length, /i/?, range between 16" and 91", and in physical coordinates betw een 0.9 and 7 
kpc. For Zq, the range i s 2" to 24" and 300 to 1100 p c. For comparison the Milky Way has ft, ^ 3 kpc l)Klvpin et alJ 
l2002j) and zq ~ 700 pc ()van der Kruit fc SearlellT98lh . The typi cal /ip/zn values for our galaxies are ~4, similar in 
value to the Milky Way and the 31 edge-on galaxies presented in lPohlen et al.l l)2000f) . 
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